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The electromagnetic radiation pressure becomes dominant in the interaction of 
the ultra-intense electromagnetic wave with a solid material, thus the wave energy 
can be transformed efficiently into the energy of ions representing the material and 
the high density ultra-short relativistic ion beam is generated. This regime can be 
seen even with present-day technology, when an exawatt laser will be built. As an 
application, we suggest the laser-driven heavy ion collider. 



I. INTRODUCTION 

The permanent interest in the problems of the interaction of relativistically strong electro- 
magnetic radiation with matter finds broad applications in the development of new concepts 
of charged particle acceleration under space and laboratory conditions. The generation of 
high energy particles, both electrons and ions, when strong electromagnetic radiation inter- 
acts with a plasma, is a well known basic phenomenon. However in the limit of extremely 
high radiation intensity it acquires new features. Under the earth conditions, an extremely 
strong electromagnetic wave can be practically realized only with laser beams. When a 
multi-petawatt laser pulse interacts with matter, conditions could be produced that were 
imagined to occur only in astrophysical objects. This opens the way for experimental studies 
of the properties of matter under these extreme conditions. With the further increase of the 
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laser power the laser based accelerators can provide the beams of the ultrarelativistic ions, 
which makes feasible their applications in the high-energy physics. 

Direct laser acceleration of protons to relativistic energies requires intensity I p = 4.6 x 
10 24 W/cm 2 x (1/im/A) 2 , corresponding to the dimensionless amplitude a = eE/m e ujc = 
m p /m e « 1836, where E, A, and u are the electric field, wavelength, and frequency of the 
electromagnetic (EM) wave, e and m e are the electron charge and mass, and m p is the 
proton mass. In a plasma, due to collective effects, protons can gain relativistic energies at 
much less intensity, about 10 21 W/cm 2 x (1/im/A) 2 , as is exemplified in the theory of the 
strongly nonlinear hybrid electron-ion wakefield induced by a short EM wave packet with 
the dimensionless amplitude a greater than (m p /m e ) 1//2 « 43 and the Coulomb explosion of 
an overdense plasma region with the size of a few microns when a relativistically strong EM 
wave sweeps all the electrons away Q]. In general, the laser-driven ion acceleration arises 
from charge separation caused by the EM wave. Various regimes have been discussed in 
the framework of this concept: the plasma thermal expansion into vacuum Q]; the Coulomb 
explosion of a strongly ionized cluster [3] ; transverse explosion of a self- focusing channel j^J ; 
ion acceleration in the strong charge separtion field caused by a quasi-static magnetic field 

a. 

In this paper we present the regime of the high density ultra-short relativistic ion beam 
generation from a thin foil by an ultra-intense EM wave. We call this regime the 'laser piston' 
(LP). In contrast to previously discussed schemes, in this regime the ion beam generation 
is highly efficient and the ion energy per nucleon is proportional to the laser pulse energy. 
Our analytical estimation conforms to the result of three-dimensional (3D) particle-in-cell 
(PIC) simulations. In comparison with the experimental experience of present-day Petawatt 
lasers, the LP regime predicts yet another advantage of the Exawatt lasers, in addition to 
possible applications depicted in (|. 



II. RADIATION PRESSURE DOMINANT REGIME 

We distinguish the following two stages of the LP operation, (i) A relativistically strong 
laser pulse irradiates a thin foil with thickness I and electron density n e . The laser pulse 
waist is sufficiently wide, so the quasi-one-dimensional geometry is in effect. Electrons are 
quickly accelerated up to v e ~ c by the transverse electric field, El, of the laser pulse and 
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they are pushed in the forward (longitudinal) direction by the force |ev e x B^/c| ps eE^, 
where Bl is the magnetic field of the laser pulse. Assume that all the electrons are displaced 
in the longitudinal direction, then the charge separation field, E\\ = 2 r Ken f X < El, between 
the electron and ion layers does not depend on the separation distance. In this longitudinal 
field the ion energy Si = (mfc 4 + (eE^ct) 2 ) 1 ^ 2 become relativistic in a time of the order of 
t ri = (rriic/eEn). We find that the ion layer can be accelerated up to relativistic energies 
during N L laser cycles under the condition E L > 1nen e l > rriiUJc/2T[eN L . Hence, to produce 
relativistic protons in one laser cycle we need an EM wave with El > 300m e ujc/e, II > 
1.2 x 10 23 W/cm 2 x (1/zm/A) 2 , and the pulse waist must be much greater than ct T i ~ A. 
(ii) The accelerated foil, which consists of the electron and ion layers, can be regarded as 
a relativistic plasma mirror co-propagating with the laser pulse. Assume that the laser 
pulse is perfectly reflected from this mirror. In the laboratory reference frame, before the 
reflection it has the energy El oc E 2 l L, after the reflection its energy becomes much lower: 
El oc E\L ps E 2 L L/A^ 2 . Here the incidence laser pulse length is equal to L, the reflected 
pulse length L is longer by factor of 47 2 , and the transverse electric field is smaller by 
factor 47 2 , where 7 = (1 — t> 2 /c 2 ) -1 / 2 is the Lorentz factor of the plasma mirror. Hence 
the plasma mirror acquires the energy (1 — I/A'j^El from the laser pulse. At this stage 
the plasma (the electrons and, hence, ions) is accelerated due to the radiation pressure. 
The radiation momentum is transfered to ions through the charge separation field and the 
'longitudinal' kinetic energy of ions is much greater than that of electrons. We note that 
the specified intensity is close to the limit where individual electrons undergo a substantial 
radiation friction effect [7]. However, when the foil is accelerated up to relativistic energy, 
this effect becomes weaker because its strength measure, 47rr e /3a 3 A, is reduced by 27 in the 
foil reference frame (r e = e 2 /m e c 2 is the classical electron radius). 

We notice a connection of the plasma layer acceleration scheme, presented above, with a 
mechanism of the ion acceleration proposed by V. I. Veksler ;8], as well as with the "snow 
plow" acceleration mechanism revealed in Ref. j^]. Formulated in the mid 1950's Veksler's 
concept of the collective acceleration of ions in an electron-ion bunch moving in a strong 
electromagnetic wave had a great influence both upon particle accelerator technology and 
plasma physics. Up to now its direct realization was considered at moderate driving EM 
radiation intensities, when transverse instabilities impede the acceleration ^(|. As we show 
below using 3D PIC simulations, in our scheme the transverse instabilities are supressed or 
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retarded due to the following, (i) The plasma layers become relativistic quickly, during one 
or more laser wave periods, in the first stage of the acceleration. Due to relativistic effects 
the transverse instabilities grow in the laboratory frame 7 times slower than in the plasma 
reference frame, (ii) The radiation pressure causes a stretching of the plasma mirror in the 
transverse direction, so the transverse instabilities can be retarded similarly to the slowing 
down of the Jeans instability in the theory of the expanding early universe 



III. PIC SIMULATION 



In order to examine the present scheme in three-dimensional geometry, whose effects 
may play a crucial role in the dynamics and stability of the plasma layer under the action 
of a relativistically strong laser pulse, we carried out 3D PIC simulations with the code 
REMP based on the "density decomposition" scheme Q). In the simulations the laser 
pulse is linearly polarized along the z-axis; it propagates along the x-axis. Its dimensionless 
amplitude is a = 316 corresponding to the peak intensity / = 1.37 x 10 23 W/cm 2 x (1 /im/A) 2 . 
The laser pulse is almost Gaussian with FWHM size 8A x 25A x 25A, it has a sharp front 
starting from a = 100, its energy is El = 10 kJ x (1 /im/A) 2 . The target is a 1A thick plasma 
slab with density n e = 5.5 x 10 22 cm~ 3 x (1/mi/A) 2 , which corresponds to the Langmuir 
frequency u pe = 7u>. For A ~ 1 /im the laser pulse electric field is strong enough to strip 
even high-Z atoms in much shorter time than the laser wave period, thus we assume that the 
plasma is fully ionized. The ions and electrons have the same absolute charge and their mass 
ratio is rrii/m e = 1836. The simulation box size is 100A x 72A x 72A corresponding to the 
grid size 2500 x 1800 x 1800, so the mesh size is 0.04A. The total number of quasi-particles 
is 4.37 x 10 9 . The boundary conditions are periodic along the y- and z-axis and absorbing 
along the x-axis for both the EM radiation and the quasi-particles. Simulation results are 
shown in figures where the space and time units are the laser wavelength A and period 
2tv/uj. 

Figure ^ shows the ion density and x-component of the EM energy flux density (the 
Poynting vector). We see that the region of the foil corresponding to the size of the laser 
focal spot is pushed forward. Although the plasma in the foil is overcritical, it is initially 
transparent for the laser pulse due to the effect of relativistic transparency (see e. g. Ref. 
[l^). Therefore a portion of the laser pulse passes through the foil. Eventually it accelerates 
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FIG. 1: The ion density isosurface for n = 8n cr (a quarter removed to reveal the interior) and 
the x-component of the normalized Poynting vector (e/m e ujc) 2 E x B in the (x, y = 0, z)-plane at 
t = 40 x 2vr/w. 

electrons and, as a result of the charge separation, a longitudinal electric field is induced. 
This can be interpreted as a rectification of the laser light, by the analogy with a recti- 
fier in electrical engineering: the transverse oscillating electric field is transfromed into a 
longitudinal quasistatic electric field. The dimensionless amplitude of the longitudinal field 
is aii ~ 150 corresponding to En = 4.8 x 10 14 V/m x (l/j,m/X). The typical distance of 
the charge separation is comparable with the initial thickness of the foil and is much less 
than the transverse size of the region being pushed. The ion layer is accelerated by this 
longitudinal field. This stage corresponds to the first step of the LP scheme. As the foil 
moves more and more rapidly, in its proper frame the incident wavelength increases, thus 
the accelerating foil becomes less transparent with time. 

As seen in the cross-section of the Poynting vector in Fig. ^ the thickness of the red 
stripes, corresponding to half of the radiation wave length, increases from left to right 
(along the x-axis). The increase is weaker at the periphery (in the transverse direction); 
correspondingly, we see distinctive phase curves. This 'anisotropic red shift' results from 
the relativistic Doppler effect when the laser pulse is reflected from the co-propagating 
accelerating and deforming relativistic mirror. The red shift testifies that the laser pulse 
is expending its energy for the acceleration of the plasma mirror, as specified above in the 
second stage of the LP scenario. The foil is transformed into a "cocoon" where the laser 
pulse is almost confined. The accelerated ions form a nearly flat thin "plate" with high 
density, Fig. El 
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FIG. 2: The isosurface for n = 2n cr , green gas for lower density at t = 100 x 2tt/lo; the black curve 
shows the ion density along the laser pulse axis. 
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FIG. 3: The maximum ion kinetic energy versus time and the ion phase space projection (x,p x ) 
at t = 80 x 2it/lu. 

Figure El shows the ion maximum energy versus time and the ion phase space plot. The 
dependence is initially linear; at later times it scales as t 1//3 . The ion and electron energy 
spectra and transverse emittances are presented in Fig. EJ The number of ions in the plate 
is Mi = 2 x 10 12 , their energies are from 1.3 GeV to 3.2 GeV. The efficiency of the energy 
transformation from laser to ions in this range is 10%; efficiency for all ions is 30%. The ion 
bunch density is 3 x 10 21 cm -3 , its duration is 20 fs, its transverse emmittance is less than 
0.1 7T mm mrad. As the plate is quasineutral, the average longitudinal velocity of the electron 
bulk is about that of ions, v e \\ ~ Vi\\, thus the average longitudinal energy of electrons is of 
the order of 8 e \\ ~ (m e /mi)£i. Correspondingly, the energy spectrum of electrons is peaked 
at much less energy than that of ions. 

According to our 2D and 3D simulations at lower intensities, corresponding to Petawatt 
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FIG. 4: The energy spectrum (red) and transverse emittance (blue) of ions (solid) and electrons 
(dashed) located in the box 50A < x < 80A, —A < y, z < A at t = 80 x 2tt/uj. The hatched region 
containes 2.7 x 10 10 //m -2 particles per cross-section. Values correspond to A = l//m. 
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ti^Petawatt pulses, the interaction exhibits a continuous transition from regimes of 
to the LP regime as the intesity increases. 



IV. ION MAXIMUM ENERGY AND EFFICIENCY 

Here we estimate the ion maximum energy and the acceleration efficiency in the model of 
the flat foil driven by the EM radiation pressure, as described above in the LP scenario. In 
general, the radiation pressure on the foil depends on its reflectance ^j]. Assume that in the 
instantaneous reference frame, where the foil is at rest, the relative amplitudes of reflected 
and transmitted waves are p and r, respectively. Here |p| 2 + |r| 2 = 1 because of energy 
conservation. The radiation pressure is the sum of the incident, reflected and transmitted 
EM wave momentum fluxes, P = (^ 2 /4vr)(l + |p| 2 - |r| 2 ) = (£ 2 /27r)(u//u;) 2 |p(u/)| 2 . Here 
the primed values correspond to the moving reference frame, variables without a prime are 
taken in the laboratory frame. In a quasi-one-dimensional geometry, at the foil location x(t) 
the laser electric field is El = EL{t — x(t)/c). If the foil is accelerated as a single whole, in 
its reference frame the incident radiation frequency becomes smaller and smaller with time, 
thus the reflection becomes more and more efficient and |p(a/)| 2 becomes closer to unity. 
In fact, the foil reference frame is not inertial since the foil is accelerated. Hence, the EM 
wave frequency lu' decreases with time in this frame as described in Ref. j^|. Nevertheless 
we can assume that the acceleration is relatively small and thus (uj'/uj) 2 = (c — v)/(c + v), 
where v = dx/dt is the foil instantaneous velocity. 

As is well known, the EM radiation pressure is a relativistic invariant jl^ . therefore we 
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can write the foil motion equation as 

dp _ E\(t - x(t)/c) | f Ap V^gg + V 2 ~ V m 
dt 2-ktiJ lP{ n yjm^ + p 2 + p ' lJ 

where p is the momentum of ions representing the foil. In the simplest case, when El = const 
and \p\ 2 = 1, the solution p{t) is an algebraic function of t. For the initial condition p = po 
at t = it can be written in a compact form as p = mjc(sinh(-u) — csch(w)/4), where 
u = (l/3)arcsinh (Qt + h^/2 + 3h /2), csch(u) = 1/ sinh('u), Q = 3El/27m e lrriiC and h = 
Po/miC + (1 +p 2 l /m 2 c 2 ) 1 / 2 . The ion kinetic energy is £jkin = m^c 2 (sinh(w) + csch(w)/4 — 1). 
As t — ► oo it asymptotically tends to ^km ~ m.jC 2 (3^£t /&im e lrriic) 1 ^ . This motion is 
analogous to that of a charged particle driven by a radiation pressure Q|, but in our case 
the role of the Thomson cross-section is played by the quantity 2/n e l. 

To find an upper limit of the ion energy acquired due to the interaction with a laser 
pulse of finite duration, we must include the dependence of the laser EM field on space 
and time. Because of the foil motion, the interaction time can be much longer than the 
laser pulse duration t^. Therefore it is convenient to consider the dynamics in terms of the 
dimensionless variable 

rt-x(t)/c 

^ J oo Aixn e lmiC d ^ 
which can be interpreted as the normalized energy of the laser pulse portion that has been 
interacting with the moving foil by time t. Its maximum value is max{?/>} = El I 'MiirtiC 2 , 
where El is the laser pulse energy, M% is the number of ions in the foil. The solution of Eq. 
rewritten in terms of ip, gives the ion kinetic energy 

2 (2x^+fr -i) 2 i r* 2 

Em n =miC — — — , k=- I \p[u)\ dtp, (3) 



2(2^ + h ) ' V jo 

for the initial condition £i(0) = (m 2 c 4 + p^c 2 ) 1 ^ 2 . The upper limit of the ion kinetic energy 
and, correspondingly, the laser to ion energy transformation efficiency can be found from 
Eq. (JH1) substituting ip = max{-0}: 

max{£ ikin } = — — — — = — T , 4 

2hEl + AlirriiC 2 Ni 

where we set po = for simplicity. Here k is the reflection coefficient, taken in the co- 
moving reference frame, averaged over the foil motion path; < x < 1. We see that if 
El N i rriiC 2 /2, in this model almost all the energy of the laser pulse is transformed into 
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ion kinetic energy. Using Eq. (J3J) and the t 1//3 asymptotic dependence of the ion energy 
on time, for given simulation parameters we estimate the acceleration time and length as 
^acc ~ (2/3)(^/A/"jmiC 2 ) 2 ti = 16 ps and x acc ~ ct acc = 4.8 mm, respectively. In the 
presented simulations the acceleration time is limited by computer resources. However, the 
analytical estimation Eq. (£Q) allows us to conlude that at given simulation parameters the 
limiting ion kinetic energy is 30 GeV. Since the ion bunch is relativistic, another ultra-intense 
laser pulse, sent with proper delay, can accelerate the bunch further. 



V. CONCLUSION 



In conclusion, the LP regime of ultra-intense laser-plasma interaction can be employed in 
a laser-driven heavy ion collider. In the collision of two ion bunches the number of reactions 
with cross-section a is JV = aAf^/s, where s is the bunch sectional area. Adopting the 
presented simulation parameters, we obtain JV k 2 x 10 30 cx/cm 2 . Provided that the energy 
is high enough so that a ~ 10 -24 cm 2 , we can get about a million events in a few femtosecond 
shot. As suggested by two of co-authors (T. T. and G. M.) in Ref. one can get a short 
multi-Exawatt laser pulse with sufficient contrast ratio (10 -12 ) using the megajoule NIF 
facility and present-day technology. Then the resulting ion bunch en ergy can be over 100 
GeV per nucleon, which is suitable for the quark-gluon plasma studies [l6j . The laser piston 
regime, being one of the examples of what we call the Relativistic Engineering, can give us 
a promising and unique tool for nuclear physics research. 
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